Multiple Andreev reflections and enhanced shot noise in diffusive SNS junctions 
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We study the dc conductance and current fluctuations in diffusive voltage biased SNS junctions with 
a tunnel barrier inside the mesoscopic normal region. We find that at subgap voltages, eV < 2A/n, 
the current associated with the chain of n Andreev reflections is mapped onto the quasiparticle 
flow through a structure of n + 1 voltage biased barriers connected by diffusive conductors. As a 
result, the current-voltage characteristic of a long SNINS structure obeys Ohm's law, in spite of 
the complex multiparticle transport process. At the same time, nonequilibrium heating of subgap 
electrons produces giant shot noise with pronounced subharmonic gap structure which corresponds 
to stepwise growth of the effective transferred charge. At eV — > 0, the shot noise approaches the 
magnitude of the Johnson-Nyquist noise with the effective temperature T* — A/3, and the effective 
charge increases as (e/3)(l + 2A/eV), with the universal "one third suppression" factor. We analyse 
the role of inelastic scattering and present a criterion of strong nonequilibrium. 
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Current transport through mesoscopic resistive ele- 
ments (tunnel barriers and disordered normal conduc- 
tors) attached to superconductors is a subject of perma- 
nent interest and intensive experimental studies. The in- 
vestigations of superconducting junctions are primarily 
focused on the complex nonlinear behavior of current- 
voltage characteristics, which exhibit subharmonic gap 
structure, zero bias anomaly, etc. In this Letter, we 
will discuss an elementary mesoscopic superconducting 
structure where the current shot noise manifests anoma- 
lous transport properties while the average current shows 
perfect Ohmic behavior. 

The circuit under discussion consists of a low-transmis- 
sion tunnel junction (or point contact) connected to volt- 
age biased superconducting reservoirs via diffusive nor- 
mal leads (Fig. la), like e.g. a short-gate Josephson field 
effect transistor [Q . In a NIN structure connected to nor- 
mal reservoirs, the average current obeys Ohm's law and 
the current fluctuations show full Poissonian shot noise, 
S = 2el, if the tunnel resistance R dominates over the 
resistance of the normal leads j| . The effect of the super- 
conducting reservoirs, which has recently attracted much 
attention, is to modify the density of states and to create 
a gap Eg in the electron spectrum of the normal leads || . 
This proximity effect provides dc Josephson current flow 
and, simultaneously, blocks the single-particle tunneling 
at applied voltages eV < 2E g . At these subgap voltages, 
the current is due to multiparticle tunneling (MPT) [Q. 
The MPT regime is manifested by the stepwise decrease 
of the current with decreasing applied voltage (subhar- 
monic gap structure at eV = 2Eg/n), which provides 
exponential decay of the current B. At the same time, 
the current shot noise undergoes enhancement due to the 
growth of the elementary tunneling charge ne || . MPT 
has been extensively studied theoretically in quantum 
point contacts and in short diffusive constrictions 

m with a wide proximity gap of the order of the energy 



gap A in the superconducting reservoirs. Both the sub- 
harmonic gap structure and the enhanced shot noise have 
been observed experimentally J9|,[l0| . 

A distinctly different transport regime occurs in long 
diffusive SNS junctions with a small proximity gap of the 
order of the Thouless energy, E g ~ Eth -C A. In this 
case, the Josephson current is suppressed, and single- 
particle tunneling dominates at virtually all applied volt- 
ages. However, when the inelastic mean free path exceeds 
the distance between the superconducting reservoirs, cur- 
rent transport at subgap voltages eV < 2A is still non- 
trivial: the tunneling electrons must undergo multiple 
Andreev reflections (MAR) before they may enter the 
reservoirs [jll[ . We will show that in long SNS structures 
with opaque tunnel barriers, the current- voltage char- 
acteristic is perfectly linear and structureless, while the 
current shot noise is greatly enhanced and reveals sub- 
harmonic gap structure (kinks) at voltages eV = 2A/n 
(incoherent MAR regime). 

The origin of the linear current-voltage dependence 
and the significant deviation of tunnel shot noise from 
the Poisson law can be qualitatively explained in the 
following way. In order to overcome the energy gap at 
low voltages (eV <C 2A) the electron has to undergo a 
large number (M sa 2A/eV) of Andreev reflections, gain- 
ing the energy eV in each passage of the tunnel barrier 
(Fig. lb). Thus, MAR tunneling in real space is associ- 
ated with probability current flow along the energy axis 
through a structure of M + 1 tunnel barriers with the 
total effective resistance Rm — {M + l)R. Since only 
electrons incoming within the energy layer eV below the 
gap 2 A participate in MAR transport, the total prob- 
ability current is I p = V/Rm- However, each pair of 
consequtive Andreev reflections transfers the charge 2e 
through the junction, and the real current I is there- 
fore M + 1 times greater than the probability current: 
I = (M + l)I p = V/R. The current flow in energy space 
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generates shot noise S p which is related to the proba- 
bility current as S p = (2/3)e7 p in the limit M — > oo. 
The arguments for the 1 /3 suppression of the Poissonian 
noise in multibarrier tunnel structures are similar to the 
ones presented in Rcf. Since the noise spectral den- 
sity is given by the current-current correlation function, 
the real shot noise S is (M + l) 2 times greater than S p , 
i.e. approaches a constant value S — (4/i?)(A/3). This 
coincides with the exact result, Eqs. (pd|), ( p"2| ) below, in 
the limit eV — > 0. 
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FIG. 1. a) Diffusive mesoscopic junction with a tunnel bar- 
rier at x — (dashed line) . b) Schematic picture of incoherent 
MAR. An equilibrium electron incoming from the supercon- 
ducting reservoir consequently tunnels through M Andreev 
side bands and creates elementary probability current Q flow- 
ing along the energy axis between the equilibrium regions out- 
side the gap. Each crossing of the barrier is associated with 
the transfer of the elementary charge e, and the electric cur- 
rent is eQ(M + 1). 

For a quantitative treatment of incoherent MAR, we 
consider the diffusive kinetic equation Jl2| for the 4x4 
supermatrix Keldysh-Green's function G{x,ti 1 t2) 



iVJ = [H, G] 



G 



1, J = VGVG, 



where V is the diffusion constant. We apply Eqs. (|l|) to 
the electrons in the normal leads and match the Green's 
functions at the tunnel barrier (x — ±0) of low trans- 
parency (r <C 1) using the boundary condition Jl3| 

J(-0) = J(+0) = (i/2)[G(-0),6(+0)], j = (l/2)v F T. (2) 

Equation (0) represents the conservation law for the su- 



permatrix current J and connects it with the voltage- 
induced imbalance of elementary probability currents j 
of tunneling electrons. Since we have assumed that the 
barrier resistance R dominates over the resistance Rd of 
the leads (R 3> Rd), the voltage drop at the barrier deter- 
mines the time-dependent phase difference between the 
reservoirs by virtue of the Josephson relation. A gauge 
transformation allows us to remove the time dependence 
in the Hamiltonian in Eq. (|l|) and cancel the electric po- 
tential. Simultaneously, a periodic time dependence ap- 
pears in the boundary condition of Eq. (|^) , which implies 
that the Green's function G(x,ti,t2) consists of a set of 
harmonics G(x, E n , E m ), E n = E + neV. The problem is 
solvable in the limit of weak Josephson coupling, r« 1 
and/or Etk <C A, when off-diagonal Green's function 
harmonics can be neglected and the diagonal harmonics 
(ri = m) satisfy the static equations, similar to the case 
of zero applied voltage. 

The dc current I can be expressed, to first order in 
r, through the quasiparticle distribution function f(E) 
defined by the following representation of the Keldysh 



function: g K = g R f-fg A , f = l(l-f) + a z f z , where g 
are the retarded (advanced) Green's functions, giving 
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dEN(E)N(E+eV)[f(E)-f(E+eV)]. (3) 



-g A ), is cal- 



The density of states, N(E) = (l/4)Trcr z (3 K - 
culated for a nontransparent barrier and normalized by 
the normal-electron density of states; all values are taken 
at the interface, x = 0. Eq. (^) has the same form as the 
conventional equation of the tunnel model ||l4|| , with the 
nonequilibrium distribution function f{E, x) obeying a 
kinetic equation following from Eq. ([!]), 



dx ^ ' ^ dx 



N(E, x) 



[f(E,x)-f (E)], 



(4) 



D(E,x) = (X>/4)Tr(l - g H g A ). The inelastic scattering 
term in Eq. (Q), describing relaxation to equilibrium pop- 
ulation fo(E) = 2np{E), is written for simplicity in the 
relaxation time approximation. The boundary condition 
for the function / at the tunnel interface obtained from 
Eq. (0) reads 



-D(E,x)(df/dx)\ = Q + (E) - Q-(E), 



(5) 



(1) Q ± (E) = ±(j/2)N(E)N(E±eV)[f(E)-f(E±eV)}. (6) 



The quantities Q , which also determine the current 
in Eq. (||), can be interpreted as spectral quasiparticle 
currents, i.e. probability currents flowing upwards in en- 
ergy space (Fig. 1): the current Q + exits from energy E 
towards energy E + eV, while the current Q~ arrives at 
energy E from energy E — eV. Along this line of rea- 
soning, the boundary condition, Eq. (|5|), represents the 
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detailed balance between the spectral quasiparticle cur- 
rent and the leakage current [the term on the left hand 
side of Eq. (0)] due to either inelastic relaxation or escape 
into the reservoirs. 

Let us consider the limit of infinitely large inelastic re- 
laxation time pl| . In this case, the leakage current is spa- 
tially homogeneous according to Eq. ([I|) . Within the en- 
ergy gap, \E\ < A, the diffusion coefficient D(E,x) turns 
to zero in the superconducting reservoirs, and therefore 
the leakage current is blocked, indicating complete An- 
dreev reflection. Thus, the spectral current Q ± is con- 
served within the superconducting gap: Q + = Q~. This 
equation provides recurrence relations for the nonequi- 
librium distribution functions f(E n ) in different side 
bands associated with MAR. The boundary conditions 
are established by the requirement of equilibrium out- 
side the gap, f(E n ) = 2n F (E n ), \E n \ > A. Indeed, 
the reservoirs maintain the equilibrium at the NS bound- 
aries, f(E,±d) — 2n F (E); on the other hand, the gra- 
dient of the distribution function given by Eq. (|^) is 
small, ~ Rd/R <C 1, and may be neglected. We note 
that the latter condition is equivalent to a small ra- 
tio between the diffusion time through the normal lead, 
d 2 /V = E^ h , and the inverse tunneling rate, (Tv F / d)^ 1 , 
i.e. Tv F /d < E Th < A. 

The physical picture of MAR in diffusive SNINS sys- 
tems is illustrated in Fig. lb. The equilibrium electron- 
like quasiparticles incoming from the left electrode with 
energy — A — eV < Eq < — A create a probability current 
Ql(E ) = jn F (E )N(E )N(E + eV) across the tunnel 
junction into the subgap region. Due to low transmis- 
sivity of the barrier and fast electron diffusion through 
the normal leads, the particle undergoes many Andreev 
reflections from the superconductor before the next tun- 
neling event will occur and, therefore, the electron and 
hole states at energy E\ = Eq + eV are occupied with 
equal probability (l/2)f(E 1 ) |l(J. Thus, the popula- 
tion f(Ei) produces both the current of holes (Ei) = 
(j /2)f(Ei)N(Ei)N(E 2 ) moving upwards along the en- 
ergy axis to the next side band, and the counter current 
of electrons Q l (Ei) = (j /2)f{E 1 )N(E )N(E 1 ) down to 
the initial state, determining the net probability current 
Q(Eq) — Q" i (E ) — Q^(Ei), and so on. As a result, the 
electron tunneling in real space is associated with the 
flow of spectral current Q + (E m ) — Q~(E m ) — Q(E ) 
through M + 1 tunnel barriers connected in series by a 
number M(Eq) = [(A — Eo)/eV] of Andreev side bands 
([x] denotes the integer part of x). In this transport 
problem in energy space, there is an effective bias volt- 
age (M + l)eV drop between the reservoirs represented 
by the spectral regions outside the energy gap, \E\ > A, 
and it is equally distributed among the tunnel barriers. 
Therefore, the distribution function has a steplike form, 



m— 1 

Z m (E Q ) = N- 1 (E k )N- l (E h+1 ), Z = 0. (8) 

k=0 

The tunnel current in Eq. (^) is determined by the 
spectral current Q + . At low temperature, the equilib- 
rium spectral current at \E\ > A is exponentially small, 
and the main contribution to the total current comes 
from the nonequilibrium subgap region. Dividing it into 
pieces of length eV and taking into account spectral cur- 
rent conservation, one finds from Eqs. (Q) 

1 f~ ^ M 4- 1 
HV)=-5 dE ——(n F (E )~n F (E M+1 )). (9) 

e-Kj-A-eV AM+1 

Equation @ describes the single-particle current in a 
tunnel junction of arbitrary length. If some MAR chain 
contains a side band E n within the proximity-induced 
gap 2E g , the corresponding density of states N{E n ) is 
zero, and the spectral current associated with this chain 
is blocked. At eV < 2E g , any MAR chain has at least 
one side band within the gap, and the total single par- 
ticle current in Eq. ([)]) vanishes. In the limit of a long 
junction, the proximity gap closes and the local density 
of states becomes constant, N(E) = 1. In this case, the 
current in Eq. (Q) shows Ohmic behaviour, / = V/R, 
with the same resistance R as in the absence of super- 
conducting "mirrors" . 

Let us turn to calculation of the tunnel current shot 
noise power S(V). A general quantum equation for the 
shot noise in superconducting junctions has been derived 
in |l7j . Assuming the asymptotic limit of highly resistive 
tunnel barrier and the long-junction approximation we 
write it on the form 

S(V)= —(f(E) + f(E+eV)-f(E)f(E+eV)). (10) 

J-oo R 

Taking into account the distribution function in Eq. (Q), 
the noise power at zero temperature becomes 



-A 



dE 



I") 



At voltages eV > 2 A this formula gives conventional Pois- 
sonian noise S = 2el. At subgap voltages, the noise power 
undergoes enhancement: it shows a piecewise linear volt- 
age dependence, dS/dV = (2e/3i?)(l + 4/([2A/eV] + 2)), 
with kinks at the subharmonics of the superconducting 
gap, eV n = 2A/n (see Fig. 2). At zero voltage, the noise 
power approaches the constant value S(0) = (4/i?)(A/3), 
which corresponds to the thermal Johnson-Nyquist noise 
with the effective temperature T* = A/3. 

The enhancement of the shot noise power can be alter- 
natively interpreted as an increase of the effective charge 
q(V) = S(V)/2I with decreasing voltage, 



f(E m ) = 2 



[n F (E M +i ) - n F (E )} - 



- + n F (E ) 



(7) 
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(12) 
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In the limit eV — > the effective charge increases as 
q(V)/e w (1/3)(1 + 2A/eV). This result differs by a 
factor 1 /3 from the value expected from a straightforward 
MAR argument which assumes the shot noise to be equal 
to the Poisson noise enhanced by the factor M. We stress 
that the 1/3 factor here results from multiple traversal 
of the tunnel barrier due to incoherent MAR and has 
nothing to do with the diffusive normal leads. 

A more detailed analysis with account of inelastic scat- 
tering shows that the current shot noise is suppressed at 
low voltage when the lifetime of a quasiparticle within 
the normal leads becomes comparable to the inelastic 
relaxation time (a > 1, see Eq. (|l4|)). Generalized re- 
currences for the distribution function in long junctions 
(N(E) = 1) then take the form pf 

f{E)-f (E) = W e [f(E+eV)+f(E-eV)-2f(E)]. (13) 

The level of nonequilibrium of the subgap electrons is 
controlled by the parameter W e = {Rd/ R){EThT £ /4). 
The strong nonequilibrium state discussed above is only 
possible at W e > 1 whereas in the opposite limit, W e <C 
1, the normal leads may always be considered as reser- 
voirs, and the enhanced noise disappears. 




FIG. 2. Spectral density S(V) of current shot noise and 
effective transferred charge q(V) = S(V)/2I as functions 
of applied voltage V. In the absence of inelastic collisions 
(solid lines), the shot noise power approaches the finite value 
5(0) = (4/i?)(A/3) at eV -> 0, and the effective charge in- 
creases as q(V) — (e/3)(l + 2A/eV). The effect of inelas- 
tic scattering is represented by dashed lines for the nonequi- 
librium parameter W £ = 5. The dependence S(V) contains 
kinks at the gap subharmonics, eV = 2A/n, as shown in the 
inset. 

Numerical results for W e = 5 are presented in Fig. 2 by 
dashed curves. The rapid decrease of S(V) at low voltage 
described by the following analytical approximation, 



^/-rr\ ^/„s 3 / , a a — sinha 
S(V) = S(0)-( tanh - + . 

a \ I smh a 



a - 



(14) 



In conclusion, we have studied subgap tunnel current 
and current shot noise in diffusive SNINS structures. We 
found that in junctions with normal leads which are much 
longer than the coherence length but much shorter than 
inelastic mean free path, the strongly nonequilibrium dis- 
tribution of the subgap electrons created by MAR is man- 
ifested in the shot noise rather than in the tunnel current. 
While the tunnel current obeys Ohm's law, the current 
shot noise is significantly enhanced and shows subhar- 
monic gap structure. 
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